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We have presented a simple approach to quantum theory of Brownian motion and barrier crossing 
dynamics. Based on an initial coherent state representation of bath oscillators and an equilibrium 
canonical distribution of quantum mechanical mean values of their co-ordinates and momenta we 
have derived a c-number generalized quantum Langevin equation. The approach allows us to im- 
plement the method of classical non-Markovian Brownian motion to realize an exact generalized 
non-Markovian quantum Kramers' equation. The equation is valid for arbitrary temperature and 
" l". friction. We have solved this equation in the spatial diffusion- limited regime to derive quantum 

' Kramers' rate of barrier crossing and analyze its variation as a function of temperature and friction. 

I While almost all the earlier theories rest on quasi-probability distribution functions (like Wigner 

^Sj ■ function) and path integral methods, the present work is based on true probability distribution func- 

' tions and is independent of path integral techniques. The theory is a natural extension of the classical 

Q , theory to quantum domain and provides a unified description of thermal activated processes and 

■ tunneling. 

; PACS number(s): 05.40.-a, 02.50.-r 
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lO . I. INTRODUCTION 

O ; 

' Ever since Kramers reported his seminal work on the Brownian motion in phase space, the theory of noise- 

■ induced escape from metastable states has become a central issue in several areas of physical and chemical sciences. 
I ' The escape is governed by Brownian motion in addition to the characteristic dynamical motion of the system in 

■ presence of a potential V{X), Brownian motion being due to the thermal forces, which in turn are associated with the 
: dissipation through fluctuation-dissipation relation at a finite temperature T. The problem and many of its variants 

have been addressed by a large number of workers over the last several decades at various levels of description 
^ I ^ I and have been extended to semi-classical and quantum domains p|-p^. A major impetus in the development of 

■ quantum theory of dissipative processes was the discovery of laser in sixties followed by significant advancement in 
^ the field of nonlinear and quantum optics in seventies and eighties when the extensive applications of nonequilibrium 

^ . quantum statistical methods were made |15 1^. Various nonlinear optical processes and phenomena were described 
^ with the help of operator Langevin equations, density operator methods and the associated quasiclassical distribution 
^ functions. However these dynamical semigroup methods of quantum optics could not gain much ground in the 
theory of activated rate processes due to the fact that these are primarily based on system-reservoir weak coupling 
^ . , and Markov approximations [ p5[ , which are often too drastic in the situations pertaining to chemical dynamics and 
' condensed matter physics. Subsequent to these developments quantum Brownian motion [p|,^|J^,p^ emerged as a 
subject of renewed interest in early eighties when the problem of dissipative quantum tunneling was addressed by 
Leggett and others and almost simultaneously quantum Kramers' problem and some allied issues attracted serious 
attention of a number of physical chemists. We refer to Ppp^ for an overview. 

The aforesaid development of the quantum theory of Brownian motion essentially rests on the method of functional 



X 



integrals |20|. This is based on the calculation of the partition function for the Hamiltonian of the system coupled 
to its environment and a non-canonical quantization procedure. The classical theories on the other hand rely on the 
partial differential equations describing the evolution of the probability distribution functions of the system both in 
the Markovian and non-Markovian regions ^ . The methods of the classical and the quantum theories are 

thus widely different in their approaches. The question is whether there is any natural extension of the classical 
theory to the quantum domain. For example, one might ask what is the quantum analogue of the classical Kramers' 
equation or Smoluchowskii's equation ? Is it possible to generalize Kramers' method of treatment of barrier crossing 
dynamics within a quantum mechanical framework so that the usual classical thermal activated process and quantum 
tunneling can be described within an unified scheme? 
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We intend to address these issues in the present paper. Specifically our object is two fold: 

1) Our primary aim here is to develop a quantum analogue of classical non-Markovian Kramers' equation, which 
describes quantum Brownian motion of a particle in a force field at arbitrary temperature and coupling in terms of true 
probability distribution function rather than quasi-probability function Pq , ^ . The generalized quantum Kramers' 
equation (GQKE) reduces to its classical counterpart in the limit fi ^ both in Markovian and non-Markovian 
descriptions. The probability distribution functions remain well-behaved in the full quantum limit. 

2) While the existing methods of calculation of quantum Kramers' rate are based on path integral techniques, we 
solve GQKE for barrier crossing dynamics as a boundary value problem taking care of the full quantum nature of the 
system and the bath. The generalized rate in the spatial-diffusion-limited regime reduces to Kramers-Grote-Hynes' 
rate in the classical limit and to pure tunneling rate in the quantum limit at zero temperature. To the best of our 
knowledge the implementation of a differential equation based approach like the present one has not been tried till 
date for a full quantum mechanical calculation of the rate. 

Based on an initial coherent state representation of the bath oscillators and an equilibrium canonical distribu- 
tion of the quantum mechanical mean values of their co-ordinates and momenta, we derive a generalized quantum 
Langevin equation in c-numbers. The bath imparts classical looking c-number quantum noise which satisfies standard 
fluctuation-dissipation relation. The simplicity of the treatment lies in the fact that the quantum Langevin equation 
is amenable to a theoretical analysis in terms of the well known classical non-Markovian theory of Brownian motion 
p^ , p5t so that the GQKE assumes the form of its classical counterpart. In what follows we show that GQKE can be 
solved in the spirit of Kramers' method to calculate the quantum escape rate. This recasting of the quantum problem 
into a classical form allows us to realize the various limits of the theory on a general footing. 

The organization of the paper is as follows. We introduce the system-reservoir model and an ensemble averaging 
procedure to derive the generalized quantum Langevin equation in a c-number form in the next section. In section 
III we derive the corresponding GQKE followed by a calculation of quantum rate of escape in section IV. The key 
result is illustrated in section V by assuming a Lorentzian density distribution of bath oscillators and a specific cubic 
potential. The paper is concluded in section VI. 

II. THE GENERALIZED QUANTUM LANGEVIN EQUATION IN C-NUMBERS 

To start with we consider the standard system- heat bath model of Zwanzig form pO| . The Hamiltonian is given by 



(1) 



Here X and P are the co-ordinate and momentum operators of the Brownian particle (the system) of unit mass and 
the set {qj,pj} is the set of co-ordinate and momentum operators for the heat bath particles. The mass of the j-th 
particle is unity and Kj is the spring constant of the spring connecting it to the Brownian particle. The potential 
V{X) is due to the external force field for the Brownian particle. The co-ordinate and momentum operators follow 
the usual commutation relation 

[X,P]=ih and [qj.pj] ^ ihSij . (2) 

Eliminating the bath degrees of freedom in the usual way |p^ , ^ we obtain the operator Langevin equation for the 
particle, 

X+ f dt'l3{t-t')X{t') + V'{X)=F{t) , (3) 







where the noise operator F[t) and memory kernel /3(t) are given by 

F{t) = ^[{gj(0) - i(0)}Kj costjji -t- Pj(0)k| sinwjt] (4) 

j 

and 

/3(t) = ^ Kj cos ujjt (5) 
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where a;| = Kj and the initial variations of the heat bath variables qj{0) and Pj{0) occur in the force term F{t). The 
relevant quantum statistical average are well known 

iF)QS = (6) 

and 



1 r 
2 



1 \ ^ ^ , , llUJ 



{F{t')F{t))Q, + {F{t)F{t'))Qs =-J2K,hcj,{coth^)cosu;,{t~t') (7) 



Here {■ ■ ■)qs refers to quantum statistical average on the bath degrees of freedom. To arrive at the above relations 
one assumes that the bath oscillators are canonically distributed with respect to the bath Hamiltonian at t = 
so that for any operator O the average is 



Tr O exp -Htath/hT 

iO)Qs = (8) 



Tr exp (^-Hbath/hT 



where Hbath — [{p'j /2) + {l/2)Kj{qj — i}^] By trace we mean carrying out quantum statistical average with number 
states of the bath oscillators multiplied by the arbitrary state of the particle. Eq.(|^) is the celebrated fluctuation- 
dissipation relation (FDR). 

Eq.(^ is an exact quantum Langevin equation in operator form which is a standard textbook material |]l5| , ^ . 
Our aim here is to replace it by an equivalent quantum generalized Langevin equation (QGLE) in c-numbers. It 
is important to mention here that again this is not a new problem ||l5| , [l6| so long as one is restricted to standard 
quasi-classical methods of Wigner functions and the like. In general, however, one is confronted with serious trouble 
of negativity or singularity of these quasi-probability distribution functions in the full quantum domain. To address 
the problem of quantum non-Markovian dynamics in terms of a true probabilistic description we, however, follow a 
different procedure. Our approach here is to split up the quantum statistical averaging procedure in two distinct 
steps. We first carry out the quantum mechanical average of Eq.(|^) 

(1) + (V'{X)) + f I3{t t'){X{t'))dt' = (F(<)> (9) 







where the averaging is taken over the initial product separable quantum states of the particle and the bath oscillators 
at t = 0, |(/)){|Q;i)|a2) • • • |aAr)}. Here \4>) denotes any arbitrary initial state of the particle and \ai) corresponds to the 
initial coherent state of the i-th bath oscillator. \ai) is given by \ai) = exp(— |Q!ip/2) X]^=o('^rVv^)l'^*)' being ex- 
pressed in terms of the variables of the coordinate and momentum of the i—th oscillator {qi (0)) — yjh/[2uji) (a^ -I- a*) 
and {pj{0)) = iyjTuoi/2{ai — a*), respectively. It is important to note that {F{t)) of Eq.(||) is a classical-like noise 
term which, in general, is a non zero number because of the averaging procedure over the coordinate and momentum 
operators of the bath oscillators with respect to initial coherent state and is given by 

{F{t)) = [{(9j(0)) - (i(0))}«, coswjt + {pjiO))Ky'' sinw.t] = /(t), say . (10) 

j 

We now turn to the second averaging. To realize f(t) an effective c- number noise we now demand that it must satisfy 
and 

if it) fit')), = i^«:,?ic.,(coth^)cosa;,(t-t') (11) 

That is, f(t) is zero centered and satisfies quantum fluctuation-dissipation relation. This may be achieved if and 
only if one introduces the following canonical distribution of quantum mechanical mean values of the bath oscillators 

att^O, 



Pj = exp 



hHfe(o))-(^(o))F + fe-(o))' 

2hu}j{nj + i) 



(12) 
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so that for any quantum mechanical mean value Oj{{pj{0)) , {(<7j(0)) — (£(0))}), the statistical average is 



ms = J O,((p,(0)), {(g,(0)) - (x(0))}) P, d{p,m d{{q,iO)) - (£(0))} (13) 

Here fij indicates the average thermal photon number of the j-th oscillator at temperature T as defined by fij = 
l/[exp(?iWj/fcBr)- 1]. 

To proceed further we now add the force term V'{{X)) on both sides of Eq.(^ and rearrange it to obtain formally 

(1) + V'{{X)) + f Pit - t'){±{t'))dt' = fit) + Q{t) (14) 

"'0 

where 

Q{t) = V'{{X))-{V'{X)) (15) 
represents the quantum mechanical dispersion of the force operator V'{X) due to the system degree of freedom. Since 



Q{t) is a quantum fluctuation term Eq.(14) offers a simple interpretation. This implies that the classical looking 
generalized quantum Langevin equation is governed by a c -number quantum noise f{t) which originates from the 
quantum mechanical heat bath characterized by the properties (^) and a quantum fluctuation term Q{t) due to the 
quantum nature of the system characteristic of the nonlinearity of the potential. Q{t) can be calculated order by 
order. In Appendix-A we show how Q{t) can be calculated in the lowest order. 



Summarizing the above discussions we point out that it is possible to formulate a QGLE (14) of the quantum 
mechanical mean value of the coordinate of the Brownian particle in a field of potential V{X), provided a classical- like 
noise term f{t) due to thermal bath satisfies (|Tl|) where the ensemble average has to be carried out with distribution 
(|l2|). To realize f{t) as a noise term we have split up the quantum statistical averaging (• • •)gg into a quantum 
mechanical mean (• • •) by the explicit use of an initial coherent state representation of the bath oscillators and then a 



statistical averaging {■ ■ ■) g oi the quantum mechanical mean values with ([121). It is easy to note that the distribution of 



the quantum mechanical mean values of the bath oscillators (12) reduces to classical Maxwell-Boltzmann distribution 
in the thermal limit huj <^ k^T, i.e. exp[— (a;|{((jj(0)) — (i(0))}^ + {pj{0))'^)/2kBT]. Secondly, the vacuum term in 



the distribution (12) prevents the distribution function from being singular at T = 0. In other words the width of 



distribution remains finite even at absolute zero, which is a simple consequence of uncertainty principle. 



III. THE GENERALIZED QUANTUM KRAMERS' EQUATION 



It is now convenient to rewrite the generalized Langevin equation ( [l4[ ) of the Brownian particle in presence of an 
external force field in the form 



x + V'{x)+ / I3{t~t')x{t')dt' ^ }{t) + Q{t) 



(16) 



where we let {X) — x for a simple notational change. f3(t) is the dissipative kernel and f{t) is the zero-centered 
stationary noise due to the reservoir where 



(/(*)>. = 0, {f{t)fit'))g=ci\t-t'\)=cir) 



(17) 



Here c(t) is the correlation function which in the equilibrium state is connected to memory kernel P{t) through FDR 
of the form 



c{t — t') = — I dijj h{lo) p{uj)fiijj 



coth 



TlLO 



2khT 



cos u!{t — t') 



(18) 



Eq.(|18D is the continuum version of Eq.(g_lJ). p{u}) is the density of modes of the reservoir oscillators. In the continuum 
version j3{t) is given by 

(19) 



j3{t — t')= I duo k{ijj)p{ijj) cos bj{t — t') 
Jo 

In the high temperature limit (fc^T >> hu) one recovers the well-known classical FDR through 
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c{t - t') = hT^(t - t') (20) 

We now proceed to the solution of Eq. . One of the essential step in this direction is to linearize the potential 
V{x) in the left hand side of Eq.([l^) around the bottom of the well a.t x = xq so that V{x) = V (xq) + {1 / 2)u!q{x — xq)'^ . 
u!q refers to the second derivative of the potential V{x) evaluated at a; = xq. This together with a Laplace transform 
of Eq. ( [l6| ) leads us to the following general solution (we take xq — for the present section) 

X{t) = {X{t))s + f M„{t - T)f{T)dT (21) 

Jo 

where 

{x{t)), = v{0)Moit) + x{0)xAt) + Go{t) (22) 

and 

Goit)= f Mo{t - T)Qo{T)dT , (23) 
Jo 

xAt) = 1-^0 / Mo{T)dT (24) 
Jo 

with a;(0) and w(0) {— i(0)) being the initial quantum mechanical mean values of the co-ordinate and velocity of the 
particle, respectively. Mo{t) is the inverse form of the Laplace transform of 

+ sPis) + w5 

with 

POO 

f3{s) = / I3{t)e-''^dt (26) 



is the Laplace transform of the dissipative kernel /3(t). The subscript in Qo,Mo and Go signifies that the corre- 
sponding dynamical quantities are to be calculated around x ~ xq. The time derivative of Eq. ( pT| ) gives 

v{t) = {v{t))s + f mo{t - T)f{T)dT (27) 
Jo 

where 

{v(t))s = v(Q)mo{t) - x{0)LulMo(t) + go{t) (28) 

with 

moit) = ^Moit) and 50 = ^Go(i) (29) 
dt dt 

It is not difficult to check that Afo(t) and mQ{t) are the two relaxation functions; mp measures how the system with a 
quantum mechanical mean velocity forgets its initial value while Mq (t) concerns the relaxation of quantum mechanical 
mean displacement. 

Now using the symmetry properties of the correlation function (/(i)/(i'))s [— c{t — t') — c{t' — t)] and the solution 
for x{t) and v{t) from (|2l|) and (|^ we obtain the following expressions for the variances 

aUt)^{[xit)-{xit)),f), 

nil 

Mo{ti)dh / Mo{t2)c{h - t2)dt2 (30a) 
Jo 

mo{ti)dti / TOo(*2)c(ii - ^2)^*2 (30b) 
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aUt) = {[x{t) - {x{t))s] Ht) - {v{t)),]), = -&,,(t) 

t i-ti 
Mo{h)dh / mn{t2)c{ti - h)dt2 

"'0 



(30c) 



The expressions for variances are general and valid for arbitrary temperature and friction and includes quantum 
effects. To recover classical hmit of the variances al^{t),aly{t),a'^y{t) one has to use (|o|) instead of (|l8|) in 30(a-c). 
It must also be emphasized that 30(a-c) are the expressions for statistical variances of the quantum mechanical mean 
values x{t) and v{t) with distribution (p2|). These arc not be confused with standard quantum mechanical variances 
which are connected through uncertainty relation. 

Having obtained the expressions for statistical averages and variances we are now in a position to write down 
the quantum Kramers' equation which is a Fokker-Planck description for the evolution of true probability density 
function p(x, v, t) of the quantum mechanical mean values of co-ordinate and momentum of the particle. To this end 
it is necessary to consider the statistical distribution of noise f(t) which we assume here to be Gaussian. For Gaussian 
noise processes we define 31, Sj] the joint characteristic function p[p^p^t) in terms of the standard mean values 
and variances, 



P, t) = exp 



ip{x{t))s + ip{v{t))s - -^{al^p^ + 2al^pp + a^^p"^} 



(31) 



Using standard procedure [p4|j2^ , |3l| , |3^ we write down the quantum Kramers' equation obeyed by the joint proba- 
bility distribution p(x, v, t) which is the inverse Fourier transform of the characteristic function; 



dp{x, V, t) 
di 



+ V\x)-^ + {no{t) - No{t))^^p{x,v,t) 



d_ 

' dx 



d_ 

' dv 



92 



92 1 



(32) 



where 



ffo(t) = Go(i) 



(33a) 



(33b) 



mo(t) 



l-wl f MoiT)dT 

Jo 



Ml{t) 



(33c) 
(33d) 



No{t) 



Yo{t) [ 



go{t)mo^{l - uj'^ I Ma{T)dT) + m^Go)(0 



(33e) 



no{t) = Mo{t)-[GQ{t)mo{t)] 



(33f) 



00 (t) = (^tiiXit) + 7o(i)<„(i) + T^Kvit) 



(33g) 



(33h) 



V{x) is the renormalized potential linearised at the bottom of the well at x = 0, the frequency being (Do(i) as 
jiven by (33d). The above Kramers' equation (^2|) is the quantum mechanical version of the classical non-Markovian 
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Kramers' equation and is valid for arbitrary temperature and friction. The quantum effects a ppea r in the description 
in two different ways. First, because of the expUcit Q-dependence go{t) [see Eqs.(23) and (|33a| )] , ^oit) and No(t) 
manifestly include the effect of quantum dispersion of the system through the nonlinearity of the potential. Second, 
the quantum diffusion coefhcients (f>o{t) and ipoit) are due to the quantum mechanical heat reservoir. In the classical 
limit ga{t), fta{t), No{t) vanishes while (j>(){t) and i/^p (t ) red uce to the forms which can be obtained by using the classical 
fluctuation-dissipation relation ( po|) in (30) and (|33 j , |33h|) . In the classical limit k^T » huj Eq.(|3^) therefore reduces 
exactly to non-Markovian Kramers' equation derived earlier by Adelman and Mazo in late seventies p^ , |25t . 

It is important to emphasize that Eq.(^) retains its full validity in the quantum limit when T Q. It is also 
apparent that Go in Eq. ( |2^ ) (therefore in go{t)) demonstrates a direct convolution of the relaxation function Mo{t) 
with quantum dispersion Qo{t). This is a clear signature of the interplay of dissipation with nonlinearity of the 
potential within a quantum description. 

The decisive advantage of the present approach is noteworthy. We have mapped the operator generalized Langevin 
equation (|^) into a generalized equation for c- numbers (|l6|) and a corresponding non-Markovian Kramers' equation(p2[). 
The present approach thus bypasses the earlier methods of quasi-probability distribution functions employed widely 
in quantum optics over the decades in a number of ways. First unlike the quasi-probability distribution function, 
the probability distribution function p{x, v, t) is valid for non-Markovian processes. Second, while the corresponding 
characteristic functions for probability distribution functions are operators, we make use of the classical characteristic 
functions. Third, as pointed out earlier the quasi-distribution functions often become negative or singular in the 
strong quantum domain and pose serious problems ||l5| , |3^ . The present approach is free from such shortcomings 
since p{x,v,t) is a true probability distribution function rather than a quasi-probability function [ p8|p9| ]. Fourth, the 
generalized quantum Kramers' equation derived here is valid for arbitrary temperature and friction. 

Regarding quantum generalized Kramers' equation we further note that although bounded the time dependent 
bmctions 70 (i) , (/>o (i) , "00(0 may not always provide the long time limits. This is well-known in classical theories 
p^ , p5| . These play an important role in the calculation of non-Markovian Kramers' rate. Therefore, in general, one 
has to work out the frequency a)o(t) and friction 70 (t) functions for the analytically tractable models. In Sec V we 
shall consider one such explicit example. 

We now consider the stationary distribution of the particles near x = which can be expressed as a solution of 
Eq.(ll) fordp%/dt^O 



d 



d 



d 



■^"5 50(00)- a)o(oo)x- (f^o(oo) - Na{oo)) 



dx 



dx 



dv 



{x,v) 



d_ 

' dv 
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7"o (00)— 1; + (/.o (00)^ + 1^0(00)^-^1 p,, (a;, w) = 0. 



(34) 



where the drift and diffusion coefficients of Eq.(34) assume their asymptotic values. 
It may be checked immediately that the stationary solution of Eq. (^) is given by 









|cxp 


[ 2Do \ 


X cxp 



V{x) + xjQQ -No + 7ogo) 



(35) 



where Dq = (/)o(oo)/7o(oo); ipo, (f)o,D,o, Nq and go are the values of the corresponding quantities in the long time 
limit. Z is the normalization constant. Here V{x) is the renormalized linear potential with a renormalization in its 
frequency. 

Eq. ( pq ) is the quantum steady state distribution. It may be checked easily that in the classical Markovian limit 
the ratio Dq goes over to kbT while "00 vanishes along with goj^lo, Nq reducing (35) to the form of Maxwell-Boltzmann 
distribution function. In what follows in the next section we shall make use of the quantum distribution (^5|) as a 
boundary condition for calculation of Kramers' rate. 



IV. THE QUANTUM KRAMERS' RATE 

We now turn to the problem of barrier crossing dynamics. In Kramers' approach the particle coordinate x (which 
in our case it is the quantum mechanical mean position) corresponds to the reaction co-ordinate and its values at the 
minima of V{x) denotes the reactant and the product states separated by a finite barrier, the top being a metastable 
state representing the transition state. 

Linearizing the motion around the barrier top at x = Xb the Langevin equation can be written down as 
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X - ujI{x - Xfc) + / /3(i - t')x{t')dt' = f{t) + Qbit) 



(36) 



where the barrier frequency lu^ is defined by V{x) = V{xb) — (l/2)uj'^{x — Xb)^. Also the quantum dispersion Qb 
has to be calculated at the barrier top. Correspondingly the motion of the quantum particle is governed by the 
Fokker-Planck equation 



(37) 



where the suffix '6' indicates that all the coefficients are to be calculated at the barrier top using the general definition 
of the last section where 



Mb{s) 



1 



(38) 



which is the Laplace transform of Mb{t) and 



Xl{t) = l+cjt / Mbit')dt' 
Jo 



(39a) 



Furthermore we have 



rrib = Mb and gb{t) = Gf,(i) 



(39b) 



^b{t) = --[lnYbit)] 



Yb[t) . ^ 



1 + Mb{T)di 



(40a) 
(40b) 



^bit) = ^ [-Mb{t)mb{t)+ml{t)] 



(40c) 



Nbit) 



Ybit) 



gbit)mb^il+ujl / Mb{T)dT) + mlGb){t) 



(40d) 



nbit) = Mbit)- [Gbit)mbit)] 
dt 



(40e) 



bbit) = ulit)al^it)+%it)alit) + -alit) 



(40f) 



Mt) = ^lait) + Ibit^lvit) + ^lit)^lAt) - 



(40g) 



In the spirit of classical Kramers' ansatz we now demand a solution of Eq. (37) at the stationary limit of the 
type. 



Pstix,v) = paix,v)Cix,v) 



(41) 



with 
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Po{x,v) = exp 



{v - gbf V{x) + x{nb -Nb + %gb) 



2Dh 



(42) 



where Db = (j)b{oo)/"fb{oo)'. ipb, gb, Nb,% arc the long time hmits of the corresponding time dependent quantities 
specific for the barrier region. The exponential factor in (41) is not the Boltzmann distribution but pertains to the 
dynamics at the barrier top a.t x = Xb- Due to the presence of 56, 0(„ Nb one may easily comprehend the signature of 
quantum nature of the system while (pb,'4'b carries the effect of quantum noise due to heat bath. The distribution po 
remains finite even at absolute zero. 

Now inserting (41), in (37) in the steady state we obtain 



ov' oxov 



Db 



,i>b + Db 



X — Xb 



{Vtb - Nb) - ^bgb) 



+ lb{v- Qb) 



dC 

dv 







(43) 



We then set 



where 



u = a{x + ab) + v - Qb 
'D.b- Nb + ^Qb + XbUJb 



and with the help of transformation (44) Eq.(43) is reduced to the following form: 



((/)fe + ail)b) 



^;^fb{x + ab) + {7. + a(l + ^)}{v - gb) 



du 



Now let 



^J^Di^ '^bix + ab) + {7fe + a{l + - gb) = Xu 



where A is a constant to be determined. 
From Eq.(44) and (47) we obtain 



where 



B C 
^± = -2A^V4l+A 



A=l + j^, B = % and C = 

Db ipb + Db 



Db -2 



By virtue of the relation (47) Eq.(46) becomes 



+ Au-^ = , 



where 



A 



_4>b + ai)b_ 

The general solution of the homogeneous differential equation (50) is 



) = F2j^ exp( —)du + Fi , 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 
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where Fi and F2 arc the two constants of integration. 

The integral in Eq.(52) converges for \u\ — > cx) if only A is positive. The positivity of A depends on the sign of a; 
so by virtue of Eqns.(44) and (47) we find that the negative root of a, i. e., a_ guarantees the positivity of A since 
—Xa = C. To determine the value of Fi and F2 we impose the first boundary condition on ( 



({x, I)) — > as a; — > 00 for all v 



This condition yields 



Fi = F-2 



(^) 



1/2 



By insertion of (54) in (52) we obtain 

CW = F2 



\2a) 



Jo 



e 2 du 



(53) 



(54) 



(55) 



Since we are to calculate the current at the barrier top, we expand the renormalized potential V{x) around x = Xb 

1 



V{x) = V{xb)-^QUx-Xbf 



(56) 



Thus with the help of (55) and (56) Eq.(41) becomes; 



Pst{x = Xb, v) = F2exp ^- 



V{xb) + Xbj^lb - Nb + gbjb) 
Db + i^b 



with 



F{x,v) 



2 du 



We now define the steady state current j across the barrier as 



vpst{x = Xb, v)dv 

An explicit evaluation of the integral using (57) yields the expression for current j at the barrier by 

Ka?{ab + XbY 



(57) 



(58) 



(59) 



j = F2 



Db 



2TTDb 



■exp{- 



exp • 



^(l + Db) 2{l + ADb) 

V{xb) + Xbj^b - Nb + gblb) \ 

Db + i^b I 



where 



F{x = Xb,v) X exp[-^^^^^^|^]rft; 

J —00 



2Db 



(60) 



(61) 



Having obtained the stationary current at the barrier top we now determine the constant F2 in Eq.(60) in terms of 
the population of the left well around x = Xq. This may be done by matching the two appropriate reduced probability 
distributions at the bottom of the left well. 

To this end we return to Eq.(41) which describes the steady state distribution at the barrier top. With the help of 
(55) we write 
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Pst{x,v) = F2 exp 



V{x) + xj^b -Nb+ gblb) 
Db + A 
2 









X exp 


[ 2Db \ 


X 



We first note that as x — > cxd , m — > 00 the preexponential factor m Pstix, v) reduces to the form 



F,[...]=F,i^ — 

We now define a reduced distribution function in x as 

"+00 



1/2 



Pstix) = / Pst{x,v)dv 



Hence from (63) and (64) we obtain 
Pst{x) = 2ttF2 



£1 

A 



1/2 



X exp 



V{x) + xjflb -Nb + gb^b) 
Db + il^b 



(62) 



(63) 



(64) 



(65) 



Similarly we derive the reduced distribution function in the left well around x — Xq using(35) as (xq may be put 
zero without any loss of generality) 



1 



P^stixo) = — V27rA^ X exp 



V{xo) + a;o(r^o - A^o + ffo7o) 



(66) 



where we have employed the expansion of V{x) as V{x) = V{xa) + (l/2)iDQ(a; — xq)^ and Z is the normalization 
constant. 

We impose the second boundary condition that at a; = a;o the reduced distribution (65) must coincide with (66) at 
the bottom of the left well i. e., 



Pstixo) =p°t(xo) 

The above condition is used to determine F2 in terms of normalization constant Z of ( |35| ) 

y(a;o) + xni^Q - iVo + go7o) 



(67) 



1 f kV^ fD,^"^ 



Z \ 27r 



X exp 



£'0 + ^0 



exp 



V{xb) - \iol{xQ - xbf + xo{Q.b - Nb + gb^b) 

Db + ipb 



Furthermore by explicit evaluation of the normalization constant using the integral 

f + 00 /' + OC 



— 00 J ~OQ 



p1t{x, v)dxdv — 1 



(68) 



(69) 



where p'ltix, v) is given by (pq). We obtain 



exp 



[»o- A^o+go7o]^ 
2{Do + A)^l 

V{xo)+ xoi^o - No + go Jo) 



ry 1/2 , ,1/2 

Z = —D„' [Do + t/Ja) X exp 
ujo 



Da+Tpo 



(70) 



Making use of (70) in (68) we obtain from (60) the final expression for quantum Kramers' rate, based on flux-over- 
population method P,|33|, as 
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LOq 



2n\2nJ dI^' {Do + ^of' 
xlDb 



exp 



[»o - iVp + go7o 
2(Do + ?Ao)c^o 
Aa2(J7f,-iV,, + 5fc7f,)2 



xexp 



.9b 



/ vr \ 1/2 , 



(71) 



where we have used the relation (56) to obtain V{xq) = V{xi,) — {1/2)oj^{xq — Xb) ^-nd definition of activation energy 
as V{xb) - V{xo). 

A close look into the definitions 39(a), 40(d), 40(e) of fl, N and g in the exponential factors in Eq.(71) reveals that 
each of them is proportional to the anharmonic correction term of the potential, V"'{x) in the leading order so that, 
[r!o-A^o+5o7o]V^o - {V'"{x))y{V"{x)), [nt- Nb + gbjb? /Ci^t ^ {V'"{x))y{V"ix)f, whereas [fib- Nb+ gblb]/Cob ~ 
{V" (x)) / {V"{x)). The last ratio being the dominant contribution the expression (71) then can be simplified as (fi 
vanishes in the long time limit) 



1/2 



Db 



2'KDb 



xexp 



[Nb - gb%)V2E 
uJb{Db + i^b) 



ADb 



9b 



\2AJ 



exp 



Db + ipb 



(72) 



The above expression is the quantum Kramers' rate which is a direct generalization of classical non-Markovian rate 
valid for intermediate to strong damping regime and for arbitrary decaying correlation function and temperature. 
The derived rate thus includes the effect of tunneling in a natural way to modify the classical rate. 

To recover the classical non-Markovian expression from (71) one has to take into consideration (i) the system 
concerned quantum correction due to non-linearity i.e., Nb and gb must vanish; (ii) the heat bath noise related 
quantities like D and ip are to be calculated from the expressions of the variances (30) using c{t — t') in the classical 
limit, i. e., (po|). Under these two conditions (71) is reduced to the classical expression 



27r 



A'^ 



1 + A'^Dt 



1/2 



exp 



E 



Df. + rb 



(73) 



where the superscript 'c' signifies the classical limit of the quantum mechanical quantities like -D o(-Dfc) and i/'o(V'fc)- 
The above expression is identical in form to one derived for classical non-Markovian dynamics m3 14, 2q]. 



V. A SPECIFIC EXAMPLE: EXPONENTIALLY CORRELATED MEMORY KERNEL 

The structure of (3{t) given in ( |l9|) suggests that it is quite general and a further calculation requires a prior 
knowledge of the density of modes p{uj) of the heat bath oscillators. As a specific case we consider in the continuum 
limit 



k{uj)p{uj) = — - 



SO that P(t) takes the well known form of an exponentially correlated memory kernel |34| 



m 



-exp 



1^1 



(74) 



(75) 



where F is the damping constant and Tc refers to the correlation time of the noise. Once we get an explicit expression 
for P(t) and its Laplace transform /3(s) = r/(l -|- sTc), it is possible to make use of (|2^) to calculate M(s) and the 
relaxation function M(t) which for the present case are given by [pj] 



Mo(s) 



aos^ 



bo + Co 



(76) 
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with 



and 



respectively, where 



1 , 2 ^ r ujI 

flo = — , oo = tj„ 4 , Co = — 



Mo{t) = cle-P"* + 46-"°' sm{et + ao) (77) 



1 , _ a„(2g"-/)-gO^-e^ 



1 , flQ \/3 



27 6 2 ^ / ' V 27 



= -^[(4 - g")' + e]'^' and ao = tan"! ( 





ao -p 


aobo 


Co 


6 


2 ^ 


1 ao 


= tan^"'^ 



(78) 



and 



^•^ = -T08 + 27 + ^~-^ + T > ° 

for the present problem. 

Now making use of the expression (77) for M(t) and expression ( |T^ ) for correlation function c{t) in Eq.(30(a-c)) we 
calculate explicitly after a tedious long but straight forward algebra the time dependent expressions for the variances 
of the quantum mechanical mean values of position and momentum of the particle. These expressions are given by 



TT Jo 1 + uj'^T^ 2kbT 

and 

In Appendix-B we provide the explicit structures of J-x{^, t) and J-'v{(^, t). Since in the long time limit cr^^, ct^^, ct^^ 
vanish, the calculation of the quantities I?Oj ^'b: V'O: V'b essentially rest on the asymptotic values of cr'^xi^) and cr^„(i) 
evaluated at the barrier top or the bottom of the well. It must be emphasized that the quantities are non-vanishing 
at T = due to quantum fluctuation of the heat bath. 

We now turn to the calculation of the relevant asymptotic coefficients and Nb in the expression for the rate {Qq 
and f2f, vanish in the long time limit because of the relaxation function). Both of them are related to the convolution 
integral through the relations (40d) and gi, — Gb- Since Gb is defined as Mb{t — T)Qb{T)dT where Mb and Qb 
correspond to the barrier top (and Mq and Qo to the bottom of the potential well) we make use of the expressions for 
Mb and Mo as given by (77) along with those for Qb and Qo as shown in Appendix-A to obtain Go{t) and Gb{t) and 
their asymptotic values. Considering only the short time linearity of Gb(t) (since the quantum effect in Qb{t) has been 
taken into account for the lowest order in the Appendix-A) calculation of gb is quite straight forward. Furthermore 
the expression for Nb as given by (40d) can be simplified in the asymptotic limit to obtain 

_ MlGb{t) 

Both Qb and Nb involve the constants {5X'^)^^o and {SXSP + 5P5X)^^q which are assumed to be h/ {2uob) (minimum 
uncertainty state) and zero, respectively for the present calculation. 
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To analyze the associated non-Markovian nature of the dynamics at various temperatures it is necessary to go over 
to numerical simulation of stationary values of Dq, Db and ipb and A. These in turn, are primarily dependent in 
cTxxit) a^iid '^vv{t)j the other variances being vanishing in the long time limit. For the present purpose we assume the 
simplest form of the cubic potential of the type V{x) = —{l/3)Ax^ + Bx^ where the parameter set used is A = 0.5; 
_ r _2 -|i/3 

B — (3/4)v4 E ; the activation energy E — V{xb) — V{xq) ~ 10. The correlation time of the noise Tc is fixed at 

0.3. The temperature and the damping constant T are varied set to set. The quantities gb and Nb which incorporate 
quantum effects through the anharmonicity of the potential can be easily calculated numerically as outlined in the 
previous paragraph. In Fig. 1 we show the Arrhenius plot, i. e, the variation of Infc vs 1/T for two different values 
of damping constant, F. It is apparent that in the high temperature regime the plot exhibits linearity, which is the 
standard Arrhenius classic result. In the low temperature regime, however, one observes a much slower variation which 
is a typical quantum behaviour. To single out this low temperature behaviour, we show in the inset of Fig. 1 a clear 
dependence of the rate - a feature observed earlier in the recent past In Fig. 2(a-c) we exhibit the variation 
of the rate fc as a function of damping constant F at several temperatures. It is apparent that at a relatively high 
temperature the rate varies inversely with the damping constant while at low temperature the rate drops at a much 
faster rate. At T = the decay is exponential in nature. The quantum rate in this situation essentially corresponds 
to zero-temperature tunneling. This result is in satisfactory agreement with that of Caldeira and Leggett ||2^. The 
present theory therefore unifies the aspects of quantum tunneling and thermal noise-induced barrier crossing on the 
same footing. 



VI. CONCLUSIONS 



In this paper we have proposed a simple approach to non-Markovian theory of quantum Brownian motion in phase 
space. Based on an initial coherent state representation of bath oscillators and a canonical equilibrium distribution of 
quantum mechanical mean values and their coordinates and momenta, we have shown that it is possible to realize a 
stochastic differential equation in c-numbers in the form of a generalized Langevin equation and the associated Fokker- 
Planck equation which can be recognized as a generalized quantum Kramers' equation. The Kramers' equation is 
then employed to derive the rate of barrier crossing which includes both tunneling and thermal induced effects on the 
same footing. The main conclusions in this study are the following: 

(i) Our ensemble averaging procedure and the QGLE are amenable to theoretical analysis in terms of the methods 
developed earlier for the treatment of classical non-Markovian theory of Brownian motion. 

(ii) The proposed Kramers' equation is an exact quantum analogue of classical generalized Kramers' equation 
derived earlier in late seventies by a number of workers. Since we have dealt here with true probability Junctions 
the theory is free from the problem of singularity or negativity of quasi- classical distribution functions which is often 
encountered in Wigner equation approaches. The equation is valid for arbitrary temperature and friction. 

(iii) The realization of noise as a classical looking entity which satisfies quantum fluctuation-dissipation relation 
allows ourselves to envisage quantum Brownian motion as a quantum generalization of its classical counterpart. The 
method is based on the canonical quantization procedure and is independent of the path integral formalisms. 

(iv) The quantum Kramers' rate (Eq.71) is valid for intermediate to strong damping regime and for arbitrary tem- 
perature and decaying noise correlations. It reduces to classical non-Markovian rate and purely vacuum fluctuation- 
induced rate or tunneling in the appropriate limits. 

(v) The theory incorporates quantum effects in two different ways. The quantum nature of the system is manifested 
through the non-linear part of the potential while the heat bath imparts the usual quantum noise. It must be 
emphasized that our general analysis takes into consideration of quantum effects of all orders . 

(vi) The theory also reveals an interesting interplay of non-linearity and dissipation in the Fokker-Planck coeffi- 
cients which include quantum corrections. The variation of the rate due to tunneling and activation with respect to 
temperature and damping has been clearly demonstrated. 

The theory presented here is a natural extension of the classical theory of Brownian motion in the sense that 
the quantum Kramers' equation is classical-looking in form but quantum mechanical in its content. Also we have 
considered only the spatial diffusion limited regime in the calculation of the rate. It is worthwhile to extend the 
approach to the energy diffusion regime and further to implement other methods of treatment of classical Brownian 
motion. 
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APPENDIX A: CALCULATION OF QUANTUM DISPERSION Q 

The quantum fluctuation Q is defined in ([l5| ) as 

Q{t)^V'{{X))-{V'{X)) (Al) 

So far as the general formulation of the theory upto Sec. IV is concerned Q is taken in full. Or in other words the 
quantum Kramers' equation ( ^2[ ) or the rate Eq.(71) incorporates quantum effects due to the system in all orders. 
However in actual calculations Q has to be estimated order by order [Q. To this end we consider the lowest order 
quantum corrections. Returning to the quantum mechanics of the system in Heisenberg picture it is convenient to 
write the operators X and P as 

X{t) = {X{t))+SX 

P{t) = {Pit)) + 5P (A2) 

{X(t)) and {P{t)} are the operators signifying quantum mechanical averages and dX and SP are quantum corrections. 
By construction {6X) and (SP) are zero and 6X and SP obey the commutation relation [SX, SP] = ih. Using (A2) in 
{V'{X)) and a Taylor expansion around {X) it is possible to express Q{t) as (keeping the lowest order non-vanishing 
term) 

Q{t)^-\v"\{X)){5X\t)) (A3) 
where {X) and {SX"^) follow a coupled set of equations as given below ; 

= (^'> 

(>) = -V'{{X)) 
j^{5X^) = {SXSP + SPSX) 

^{SXSP + 5P5X) = 2{5P'^) - 2V"{{X)){5X^) 



^{SP^) = -V"{{X)){SXSP + SPSX) (A4) 
at 

The above set of equations can be derived from Heisenberg's equation of motion. In the Kramers' problem one 
is primarily concerned with the local dynamics of the system around the barrier top at {X)(= x) — Xb or at the 
bottom at {X){= x) = xq. The solution depends critically on the nature of curvature of the potential, i.e., V"({X)). 
Considering the local nature of the dynamics we have expressed V"{x — xj,) = —i^l and V"{x — Xq) = cJq and 
consequently the equations for the quantum corrections can be solved independently of the first two of Eqs.(A4). A 
solution of {5X'^) around x ~ xq which is an elliptic fixed point is therefore given by 



lr/r-2v (^^ )t=oi {5X5P + 5P5X) ^ 2/, 



{5X\t))^^^^ = -m'),^, - ^^jr^] cos{2LOot) + ^ — ^ sin(2c^oi) + ^ (A5) 
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where /q is an integration constant and is given by Jq = {SP"^) + ujQ{SX'^)^^f^. Therefore the quantum dispersion at 
the bottom of the weU is given by, 

Qoit) = -lv"'{xo){SX^{t))^^^^ (A6) 

where {^X'^{t)) ^^^^ is governed by Eq.(A5). Similarly we calculate quantum fluctuation Qb{t) near the top of the 
potential barrier at 

Qb{t)^~lv'"ix,){SX\t))^^^^ (A7) 
where {SX'^{t)) in a solution of the last three equation of (A4) around x ^ Xb and is given by 



(^^' (0) ... = {6X') cosh(2^fct) + iS^^P+JP^X)t=o sini,(2^,t) ( A8) 

ZUJb 



It is interesting to note the hyperbolic nature of the top of the barrier which is reflected in the exponential divergence 
of the quantum fluctuations. This point has been studied extensively in the recent literature in the context of chaos 
p6|-09t. Having evaluated the quantum dispersions Qo(t) and Qb{t) we are now in a position to calculate several 
related quantities like G{t), N{t), fl{t) and g{t). A better estimate of the quantum correction Q can be obtained from 
the solutions of the equations of higher order corrections derived earlier by Sundaram and Milonni | ]35[ | . 

APPENDIX B: CALCULATION OF Tx{uj,T) AND J^v{iJ,T) 
The expressions J^^(Lu,t) and !Fy{uj^t) are given by 

/ A/o(ti)[c?/i(ti) + c^/2(ti)]dti (Bl) 

and 

:F^{uj,t)^ / No{ti)[p\\h{ti) + qoclhiti) - S,^h{ti)]dti (B2) 







where 

Na{t) - c°/e-P°* + cO^Oe-""* sin(ei + a) - cOee"""* cos{et + a) (B3) 
Here /i, /2, and /a are given by the following expressions. 

/i = -B2je~^*^ + B2j cos cujti + Bij sin ujjti (B4) 

I2 — ^[—B4j cos{ujjti + a) + B4je^'^*^ cos{ujjti + a) cos(cjj_ii) + B^je^'^*^ cos{tJjti + a) 

X sin(a;j_ti) + B^j sm{ujjti + a) — i?6je~''*^ sm{ujjti + a) cos{ujj-ti) + B^je^'^*''^ x 

sin(wjti + a) sa\(ujj-ti) + B^j cos[ujjti — a) — Bj,je^'^^^ cos{u}jti — a) cos(cjj+<i) 

— Ssj-e^'*^ cos{ujjti — a) sm{u}j+ti) + B^j sin(a;jti — a) — B^je^"^^^ sin{u}jti — a) x 

cos(cL'j+ti) + B^jC^'^*^ sm{ijjjti — a) sin(cL'j+ti)] (B5) 

I3 = ^[B^j cos{ujjti + a) — Sgje^**^ cos{ujjti + a) cos{ujj^ti) + B^jC^'^*'^ cos{ujjti + a) 

X sin LUj^ti + B4j sin(a;jii + a) — _B4je~'*^ sm{ujjti + a) cos{ujj^ti) — 

-Bej-e'**^ sin(a;jti + a) sm{ujj^ti) + B^j cos{ujjti — a) — i^sj-e"**^ cos(ajjti — a) 

X cos(a;j+ii) + Bsj-e"'*^ cos(cjjii — a) sin(a;j+ii) + B^j sm{ujjti — a) — Bsj-e"'*^ 

X sin(a;jti ~ a) cos{ujj+ti) — B^je^''*^ sm{ujjti — a) sm{ujj+ti)] (B6) 
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where 



^3+ = ojj + a 



uij- = Ui — a 



Q 



B2j = 



Bij - 



P 



= 9 , n 



- w 
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FIG. 1. Plot of Ink vs 1/T using Eq.(72) for different values of the damping contant (a) F = 1.3, and (b) F = 1.7 [ Inset : 
Plot of A: vs T to illustrate the dependence at very low temperature. The parameters are same as in the main figure ] (units 
arbitrary). 



FIG. 2. Plot of quantum Kramers' rate k vs T using Eq.(72) for (a) T = 5.0, (b) T = 3.0 and (c) T = 0.0 (units arbitrary). 
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